Classical bifurcation in a quadrupolar NMR system 
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The Josephson Junction model is applied to the experimental implementation of classical bifurca- 
tion in a quadrupolar Nuclear Magnetic Resonance system. There are two regimes, one linear and 
one nonlinear which are implemented by the radio-frequency term and the quadrupolar term of the 
Hamiltonian of a spin system respectively Those terms provide an explanation of the symmetry 
breaking due to bifurcation. Bifurcation depends on the coexistence of both regimes at the same 
time in different proportions. The experiment is performed on a lyotropic liquid crystal sample of an 
ordered ensemble of ^^^Cs nuclei with spin / = 7/2 at room temperature. Our experimental results 
confirm that bifurcation happens independently of the spin value and of the physical system. With 
this experimental spin scenario, we confirm that a quadrupolar nuclei system could be described 
analogously to a symmetric two-mode Bose-Einstein condensate. 

PACS numbers: 05.45.Mt, 03.75.Mn, 03.67.-a 



in 
in 

o 



The Josephson Junction (JJ) model remains one of the 
key concepts for theoretical advances in the physics of 
superconductivity and superfluidity. Within the ultra- 
cold atom scenario, the description of two-mode Bose- 
Einstein condensates (BEC) by means of the JJ model 
has afforded new insights into nonlinear tunneling Q , 
owing to nonlinear ity as a source, revealing many kinds 
of phenomena, from entanglement to classical bifurca- 
tion [3|-[iq|. Lately, classical bifurcation effects have at- 
tracted the attention of researchers, as they can indicate, 
for the associated quantum systems, a signature of quan- 
tum phase transitions jlll-ll4j and, more recently, have 
been used in the study of an unstable quantum pendu- 
lum [iS^. In this context and due to the huge advances 
in experimental control, the investigation of systems such 
as two BEC traps |5|, |lO| or two vibrational degrees of 
freedom in polyatomic molecules [ig'-^S^ has become very 
active. The Hamiltonians modeling such systems are well 
described by raising and lowering operators and they may 
be rewritten in terms of the SU(2) operators and their 
commutation relations via the Schwinger pseudospin rep- 
resentation. In this representation, these Hamiltonians 
display nonlinear terms that commonly appear in the z- 
component, as for instance in the case of a symmetric 
trap of a two-mode BEC. One may write the simplest 
Hamiltonian model as H = ~ ^J^c, where x ^"^^ 
Vt represent respectively the nonlinear coupling due to 
atom-atom interaction and the linear coupling due to an 
external perturbation. By this and similar models, some 
efforts have been made to understand the evolution of 
the wave function \^ {t)) in the spin-coherent represen- 
tation, and to analyze entanglement chaos in kicked 
spin systems and bifurcation (ioj . 

The Hamiltonian in the Nuclear Magnetic Reso- 
nance (NMR) scenario, has a direct physical interpreta- 



tion and in the rotating frame corresponds to the Hamil- 
tonian of a quadrupolar nuclear system acted on by 
a radio-frequency pulse along the negative x-direction 
2l|, |22|. In this letter, we explore this equivalence of in- 
terpretation and study the feasibility of observing a sig- 
nature of bifurcation in quadrupolar spin systems, while 
also assessing the use of the JJ-model in nuclear systems. 

Nuclear spin systems with / > 1/2 are described by the 
Zeeman term, the quadrupolar term, the radio-frequency 
term and weak interactions with, among other things, 
nuclei, electrons and field fluctuations, which we refer to 
as an environment term denoted by He 
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The Zeeman term is the interaction between the spin 
nuclear magnetic moment —h'yl = —hj{Ix^ly,lz) and a 
strong constant magnetic field Bq = (0, 0, 5o) aligned in 
the z-direction. The quadrupolar term is due to the in- 
teraction of the nuclear quadrupole moment (Q) with an 
electric field gradient {Va(3^ with a, /3 = x^y^z)^ such that 
in an oriented system with axial symmetry, the following 
inequality is satisfied \Vzz \ ^ IV^^^I ^ \^yy\^ which allows 
us to express the term as ^j^2^^-i) i^^'z ~ -'-^)- "^^^ radio- 
frequency (RF) term represents the interaction between 
the spin nuclear magnetic moment and an external time- 
dependent magnetic field, which is perpendicular to the 
strong constant magnetic field Bq; this term is written 
as -{-hjBi {Ix cos {ojRFt -\- (j)) -\-Iy sin {ojRFt + </>)), where 
the phase (j) defines its direction on the xy-plane. In a 
rotating frame, the NMR Hamiltonian is given by 



-\-JlWl {Ix COS (j) + 
3eQV, 



ly sin ( 



where uq = 2/(2/-T)^ represents the quadrupolar cou- 
pling, uji = jBi is the RF strength, and ul 
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Larmor frequency of a nuclear species. The Larmor fre- 
quency and quadrupolar coupling satisfy the inequality 

I^lI > |^q|. 

To match the Hamiltonian 1-L and the NMR Hamil- 
tonian Hnmr^ let us choose ujrf = ujl and (/) = tt, 
such that, without loss of generality, we may drop the 
constant term —^^t^ and the environment term W^nv 
Now the NMR Hamiltonian takes the form Hnmr = 
hwilx. Next, by substituting the dimension- 



huj 



Q j2 



less parameter A 
rewritten: 



the NMR Hamiltonian can be 
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We then use this Hamiltonian to look into the classical bi- 
furcation mechanism in nuclear spin systems. The corre- 
sponding semiclassical Hamiltonian is generated by map- 
ping the quantum mechanical operators onto the com- 
plex numbers, following the usual procedure reported 
in [3I, l2o|. This amounts to letting Iz ^ z and I^; — ^ 
\/l — z'^ cosC, giving 
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W = -z^-^l-z^ cosC, (3) 

where, for our purposes, z represents the temporal mean 
z-magnetization and ^ a relative phase. In classical me- 
chanics, Eq. (|3]) describes the motion of a particle in 
a phase potential defined by V {Q = —y/l — z^ cos 
where V is shaped by cos^ and weighted by \/\ — z^ 
as sketched in Fig. [TJa). A is a tunable parameter that 
determines the dynamics of a particle in a conserved en- 
ergy configuration. From Hamilton's equations of motion 
z = —dl-L'/dC, and C = dl-L' /dz we have, 

z = -\J\-z^ sinC, (4) 

C = A^+ ^-^ cosC, (5) 
V 1 — ^ 

such that z and ^ are canonically conjugate variables. 
The fixed points of the Hamiltonian denoted by 
P = (2:0, Co), are determined by the condition i = = 0. 
These are ("o = ±^7r with n — 0,1,2,... and = 
|o, ±y^l — 1/A^|. Note that there are many interest- 
ing sets of fixed points to analyze. The trivial one, 
Po = (0,0), which corresponds to a stable fixed point, 
is located on the positive X-axis of the coordinate frame 
(see Figlljb)). Physically, in the NMR scenario, it cor- 
responds to an extreme situation of null quadrupolar 
coupling, or to a linear regime (A < 1). Still in the 
linear regime, P^^ = (0, tt) matches another stable fixed 
point, located on the negative X-axis of the coordinate 
frame (see Fig [Hb)). The non-trivial fixed points are 

P± = ( ±Y^1 — 1/A^, TT j ; here, for A > 1, it can be seen 



FIG. 1: (Color online) Sketch of the spin scenario of the bi- 
furcation process discussed in this letter, (a) Initially the 
spin state precesses around fixed point Po, such that after 
the dynamics the spin system goes to one of the points P± 
depending on the value of A. (b) Trajectories drawn on the 
spherical phase space before and after the bifurcation process 
that depends on the value of A. The typical supercritical 
pitchfork bifurcation scenario occurs; i.e., a stable fixed point 
bifurcates into two new stable fixed points while the original 
becomes unstable. 

that the stable fixed point P^^ undergoes a supercritical 
Pitchfork bifurcation, becoming unstable, and it splits 
up into two others, the P± stable fixed points. These 
are divided by a separatrix in Fig. [T]^b). In the NMR 
interpretation, this picture corresponds to the situation 
of a quadrupolar coupling stronger than the RF pulse 
intensity, or to a nonlinear regime. Our efforts were fo- 
cused on finding out how this theoretical analysis could 
be reached by using nuclear spin systems. 

The spin scenario and the experimental imple- 
mentations were performed in a lyotropic liquid 
crystal sample prepared with 42.5 wt % cesium- 
pentadecafiuoroctanoate (Cs-PFO) and 57.5 wt % 
deuterated water (D2O) [23|. Cesium nuclei (^^^Cs) are 
quadrupolar spin systems with 1 = 7/2 such that the 
dimension of the Hilbert space is d = 2/ + 1 = 8. The ex- 
periment was carried out in a Varian 500 MHz - Premium 
Shielded (11.7 T) spectrometer at room temperature (25 
°C). A liquid NMR 5mm probe was used in this experi- 
ment. The Larmor frequency and quadrupolar coupling 
are 65.598 MHz and 7.7 kHz, respectively. The 7r-pulse 
time was calibrated as 25 jas. The spin-lattice and spin- 
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spin relaxation time, for cesium nuclei, is Ti ^ 320 ms 
and T2 « 4 ms, respectively. The recycle delay time is 
di 1.8 s. 

To describe a quantum state in NMR implementa- 
tions, we use the density operator formalism representing 
the thermal equilibrium state, whose populations satisfy 
the Boltzman-Gibbs distribution. Theoretically, the den- 
sity operator is represented hy p = ^ exp [— /3Ho], where 
Ho = —hujLiz is the secular contribution of the NMR 
Hamiltonian and (3 = {kBT)~^ ^ ks being the Boltzmann 
constant T the room temperature. If the polarization 
strength is e = Pujl^/ where Z is the partition func- 
tion, and this factor has a value ~ 10~^ then the density 
operator could be expanded to a first order approxima- 
tion as p = — epo, in which po = is called the 
deviation density matrix. 

To initialize the quantum state, we transform po to 
prepare a pseudo-nuclear spin coherent state (pseudo- 
NSCS) using the protocol of Ref. [H, [25|. The pseudo- 
NSCS is denoted as so the density operator is 

p= (I -e) 1 + eAp, such that Ap = |C (^)) (C ^)|, 
for any < ^ < tt and < (/? < 27r. 

For our purpose, the angular parameters (polar and 
azimuthal if) of a pseudo-NSC S were chosen in such a way 
that the phenomenon of bifurcation appeared as sharply 
as possible and therefore these parameters were fixed at a 
pair of initial conditions |C+ (7r/4,7r)) and |C- (37r/4,7r)), 
where the ("-positive (^-negative) represents an initial 
condition at the north (south) hemisphere of a spherical 
phase space. 

To sketch the classical bifurcation we use a control 
parameter A, such that uoq is maintained at a con- 
stant strength and uji is varied from highest to low- 
est values. The different strengths of uoi are quantified 
by the calibration of 7r-pulses at various elapsed times 
=25, 30, 40, 50, 60, 100 /is or for the parameter 
A = 0.67, 0.81, 1.08, 1.35, 1.62, 2.70. Once we have cho- 
sen the strength of cji, the implemented pseudo-NSCS 
|C± (^5 ^)) is transformed by a RF pulse that depends on 
the Hamiltonian in Eq. (|2]) at various times Tk = /cAr 
and k = 0, . . . , 44 with Ar = 5/is. At each step, the 
pseudo-NSCS is tomographed by the quantum state to- 
mography (QST) procedure [H, H^. Next, the average 
value of the z-component of the spin angular momentum, 
(Iz {Tk)) = Tr {Ap(r/e) I;^}, is computed. 

In Fig. [2] we show the experimental results for z- 
magnetization and its distribution for initial conditions 
on the north (south) hemisphere of the spherical phase 
space under the linear (nonlinear) regime. On the left 
of Fig. [21(a) are shown experimental (dots) and nu- 
merical (solid line) results [27| for (I^ (r)), for the ini- 
tial condition on the north hemisphere under the lin- 
ear regime, but not too far from the nonlinear regime, 
with approximately A ^ 0.67. The beats and fast de- 
cay is a typical signature of an intermediate regime. On 
the right of Fig. [2ja) we display a histogram of the 
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FIG. 2: (Color online) Experimental results (symbols) and nu- 
merical results (continuous lines) for the dynamics of mean av- 
erage values of ^^-magnetization (left) and a z-magnetization 
distribution (right). Each value was calculated from the to- 
mographed deviation density matrix at 45 different values of 
r. Results for (a) initial condition on north hemisphere and a 
linear regime, (b) initial condition on south hemisphere and 
a nonlinear regime. 
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FIG. 3: (Color online) Experimental results (symbols) of bi- 
furcation phenomena in spin systems which are described by 
the JJ-model. Two initial conditions are denoted by dashed 
lines, which correspond to and pseudo-NSCS. The 
spin system can be driven under different regimes: a linear 
(0 < A < 1) and a non-linear (1 < A). Theoretical predic- 
tions (solid lines) are sketched by using zhy^l — 1/A, whose 
positive (negative) sign corresponds to the north (south) 
hemisphere of a spherical phase space. The cyan circle (ma- 
genta square) indicates that the data was calculated from ex- 
perimental results explained in Fig. (Fig. Elb)). The 
error bars are bounded at 10 % of the maximum value of the 
time averaged ^-magnetization, which is |=bl|. 
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experimental values of (I^ (r)) and the Gaussian distri- 
bution Gz ((I^)) = Go exp — ((I^) — zf' / crl has been 
calculated and drawn. The parameter referred to as 
the temporal mean z-magnetization, corresponds to the 
mean value of (I^) over 45 different elapsed time points. 
Similarly, the parameter cFz is the well-known standard 
deviation. The main information extracted from Gz is 
the z-value, which depends on A, and this parameter 
indicates the stage of bifurcation. 

On the left of Fig. [2l^b) we show experimental (squares) 
and numerical (solid line) results for (I^ (r)) in the non- 
linear regime, satisfying A = 2.7 > 1. The smooth beats 
are almost completely attenuated and the decay is slower 
than in the linear regime. This happens because uoi is 
weaker than ujq. On the right of Fig. [2](b), there is a his- 
togram for experimental values of (I^ (r)) and a Gaussian 
distribution ((I^))- 

In Fig. [3] we present the experimental results (sym- 
bols) and theoretical prediction (solid lines) of the tem- 
poral mean z-magnetization for initial conditions on the 
north hemisphere (dots) and south hemisphere (squares). 
The cyan circle (magenta square) indicates that the data 
was calculated from experimental results explained in 
Fig. [2fa) (Fig. Efb)). We observe that the experimental 
results match the theoretical prediction of bifurcation. 
To explain this phenomenon in a nuclear spin system, 
we need to remember that the eigenstates of the secular 
Hamiltonian are |m), with m — — /, — / + 1, . . . , / — 1, /, 
and the operator that depends on Hamiltonian ([2|) trans- 
forms each eigenstate. In the linear regime, the RF term 
imposes the dynamics of spins, maintaining the spin sys- 
tem under a superposition of the |m) eigenstates, lead- 
ing the quantum state from to and vice- versa. 
This is analogous to what happens in the tunneling phe- 
nomenon of a symmetrical trap of a two-mode BEG [H- 
H, what happens in the description of the JJ- 

model in superconductivity, called the plasma oscillation 
regime On the other hand, in the nonlinear regime, 
the nonlinear term decides the behaviour of the system. 
In this case, from the basic principles of Quantum Me- 
chanics, we know that |±m) = hw? |±m); for |+m) 
and \ —m) there is a degeneracy which drives the spin sys- 
tem to the bifurcation phenomenon, because the initial 
condition for the north (south) hemisphere is retained by 
the operator that depends on Hamiltonian ([2]) and the 
spin system is induced to precess aligned parallel (an- 
tiparallel) to an effective magnetic field which is aligned 
along the z-axis of a reference frame. Similarly, in a two- 
mode BEG, this regime corresponds to the self-trapping 
regime [I]. 

Finally, we draw our conclusions: we performed a clas- 
sical bifurcation in a nuclear spin system that is described 
and interpreted by the JJ-model. The spin scenario coin- 
cides with other experimental schemes commonly named 
as the symmetric double- well trap of two-mode BEG, 



with the possibility of extending to the asymmetric case. 
We take advantage of the main physical property of a ly- 
otropic liquid crystal sample, which is the collective order 
in the presence of a magnetic field. This inspired us to ex- 
plore the nonlinear regime to study squeezing processes, 
which are currently being developed in our laboratory. 
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